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ABSTRACT 

Polarization behaviors of optical afterglows from two-component gamma-ray burst 
jets are investigated, assuming various configurations for the two components. In most 
cases, the observed polarization is dominated by the inner narrow component for a 
long period. Interestingly, it is revealed that different assumptions about the lateral 
expansion of the jet can lead to different evolutions of the position angle of polarization. 
The observed afterglow light curve and polarization behaviors of GRB 020813 can be 
well explained by the two-component jet model. Particularly the model is able to 
explain the constancy of the observed position angle in this event, given that the line 
of sight is slightly outside the narrow component. 
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1 INTRODUCTION 

Soon after the discovery of gamma-ray burst (GRB) after- 
glows, it was expected that radiation from these GRB em- 
bers may be polarized (Loeb & Perna 1998; Gruzinov & 
Waxman 1999; Medvedev & Loeb 1999). The first direct de- 
tection of polarized afterglows from GRB 990510 by Covino 
et al. (1999) and later confirmation by Wijers et al. (1999) 
intrigued subsequent investigations of the geometry of GRB 
ejecta (Ghisellini & Lazzati 1999; Sari 1999; Gruzinov 1999). 

Polarization originates naturally from some asymme- 
tries when relativistic electrons produce non-thermal radi- 
ation. To model the observed polarizations of GRB after- 
glows, there are basically two kinds of consideration on the 
violation of symmetry with respect to the observer. The 
first class includes the models assuming ordered magnetic 
fields. The magnetic fields can be either locally ordered, 
which corresponds to the magnetic domain model (Gruzi- 
nov & Waxman 1999), or even entirely aligned within the 
ejecta which is magnetized by the central engine (Granot 
& Konigl 2003; Lyutikov, Pariev & Blandford 2003; Lazzati 
et al. 2004; Dai 2004). The second class involves beaming 
effect, where the GRB ejecta is assumed to be conically col- 
limated. Observers need to be off-axis, which is natural since 
a large viewing angle corresponds to a large chance possi- 
bility. The magnetic fields are postulated to be randomly 
but an-isotropically distributed behind the shock. The sim- 
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plest configuration is two dimensional, with the magnetic 
fields randomly distributed in the shock plane (Ghisellini 
& Lazzati 1999; Laing 1980). In further particular advise- 
ments the level of the magnetic field anisotropy has been 
parameterized (Gruzinov 1999; Sari 1999; Granot & Konigl 
2003). Previous works revealed that the simple homogeneous 
jet is confronted with abrupt variation of polarization an- 
gle by 90°, when the degree of polarization passes through 
zero. No intermediate values for the position angle are ex- 
pected in this model. In fact, the observed position angle of 
some afterglows, e.g. GRBs 990712 (Rol et al. 2000), 020405 
(Covino et al. 2003a) and 020813 (Gorosabel et al. 2004), did 
not exhibit such significant changes 1 . As indicated by sim- 
ulations of jet propagating in the collapsar, the emergent 
jet from the stellar envelope is structured (Zhang, Wooley 
& Heger 2004). Two approximations on the jet structure 
had been made previously. One is the power law structure 
(Meszaros, Rees & Wijers 1998; Dai & Gou 2001; Rossi, Laz- 



1 Significant evolution of position angles does have been observed 
in two events, i.e. GRBs 021004 and 030329 (Rol et al. 2003; 
Greiner et al. 2003). However, this does not necessarily imply the 
homogeneous jet structure since their light curves are peculiar. As 
first pointed out by Granot &: Konigl (2003), variable afterglow 
light curves are expected to be accompanied by variable polariza- 
tion light curves. Nakar & Oren (2004) fitted the light curve and 
polarization curve of GRB 021004 within the angular inhomoge- 
neous "patchy" shell model, while Bjornsson, Gudmundsson & 
Johannesson (2004) interpreted this afterglow within the context 
of the refreshed shock model. 
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zati & Rees 2002), while the other is the Gaussian profile 
(Zhang & Meszaros2002). Recently, Rossi et al. (2004) had 
detailedly investigated the polarization behavior of both a 
power law structured jet and a Gaussian jet. The advantage 
of these kinds of structured jets is that the position angle 
keeps constant with time. 

In this paper we explore the polarization of GRB af- 
terglows by assuming a different jet structure, i.e. the two- 
component jet structure. The two-component jet model was 
proposed to explain observations, such as unusual radio to 
X-ray spectrum and late time bumps in some afterglows 
(Frail et al. 2000; Berger et al. 2003; Huang et al. 2004). 
Liang & Dai (2004) found that there exist two peaks in the 
histogram of the spectral peak energy distribution derived 
from vF u in GRBs, which also tentatively implies a two- 
component structure of GRB jets. In Section 2 we describe 
the dynamics and radiation mechanism of two-component 
jets. The afterglow polarization behavior of two-component 
jets is presented in Section 3. Observed light curve and po- 
larization curve of GRB 020813 are fitted within this model 
in Section 3.3. We conclude and discuss our results in Section 
4. In the Appendix, the polarization degree of synchrotron 
radiation in an ordered magnetic field as a function of fre- 
quency is given. 



2 JET HYDRODYNAMICS AND 
SYNCHROTRON RADIATION 

For simplicity, we make several assumptions on the hydro- 
dynamic evolution of a two-component jet. First, the two- 
component jet evolves adiabatically. Second, we apply the 
uniform thin shell approximation to both components of the 
jet. Third, we assume that there is no interaction between 
the narrow and wide components. Each component evolves 
independently and therefore can be treated separately (see 
also Huang et al. 2004). Lastly, two different assumptions 
are made on the lateral expansion of the jet. One is to as- 
sume that the jet expands laterally at the speed of light, and 
the other is to postulate that the jet experiences no lateral 
expansion. The former corresponds to the case of a clear-cut 
lateral boundary between the uniform jet and the environ- 
ment, which was assumed for simplicity in previous works 
(Rhoads 1999; Sari, Piran, & Halpern 1999). However, the 
physical parameters near the actual lateral boundary should 
vary gradually (Zhang et al. 2004) . Recent hydrodynamical 
calculations showed that for smooth jet profiles the lateral 
expansion speed is much less than the sound speed in the jet 
co-moving frame even during the relativistic stage (Kumar 
& Granot 2003; Granot & Kumar 2003). So, whether the jet 
expands laterally or not is still an open question. We there- 
fore also discuss the second possibility that the jet has no 
sideways expansion. In the following, we denote the physi- 
cal parameters of the narrow and wide components by the 
subscripts "TV" and "W", respectively. The narrow (wide) 
component can be described by the isotropic kinetic energy 
-Ejv.iso (Ew,iso), the bulk Lorentz factor 7jv {fw), and the 
half-opening angle 9n (0w)- 

The hydrodynamic evolution can be evaluated in the 
same way for both components of the jet. One caveat should 
be made about the inner boundary of the wide component. 
The wide component is assumed to be a hollow cone, with 



the inner half-opening angle to be constant and the same 
as the initial aperture of the narrow component during the 
whole relativistic stage. However, the hydrodynamics of the 
wide component should be hardly affected by this fact. We 
let #o be the initial half-opening angle of an arbitrary compo- 
nent. After a short coasting phase or reverse-forward shock 
interaction epoch, the jet begins to decelerate in the inter- 
stellar medium (ISM), with its radius R, bulk Lorentz fac- 
tor 7 and half-opening angle 9j evolve as (Sari, Piran, & 
Narayan 1998) 

R = 5.8 x lO 17 ^ 4 ™- 1 / 4 ^!^)- 1 / 4 cm, (1) 

7 = ll^n-VV^^) 378 , 0j ~ 9 , (2) 

where the isotropic kinetic energy E = 10 53 i?53 erg, the ISM 
number density n is in units of cm~ 3 , td is the observer's 
time in days, and z is the redshift. During this stage, the 
hydrodynamics of the jet is the same as that of an isotropic 
fireball, which obeys the self-similar solution of Blandford- 
McKee (1976). As the jet continues to decelerate, there is 
a critical moment when the Lorentz factor 7 = 6q , which 
corresponds to the time 

tj = 0.645(1 + z)EH 3 n- 1/3 ef/\ day, (3) 

where 80 — 0.1#o,-i- In the case of no lateral expansion, the 
jet hydrodynamics after this time is unchanged and still can 
be described by equations (1) and (2). On the other hand, 
if the jet expands laterally at the speed of light, the jet will 
experience a runway sideways expansion (Rhoads 1999; Sari 
et al. 1999). The hydrodynamic quantities in this case when 
t > tj evolve as 

R = 6.18 x lO 17 £^ 3 n" 1/3 0o^i cm, (4) 

7 = * - 1 (fr 1/2 , fi«7" 1 . (5) 

Polarized afterglows have been observed usually within sev- 
eral days (Covino et al. 2002; Bjornsson 2002). In this paper, 
we focus on the stage when the jet is still relativistic. The 
moment when the jet becomes non-relativistic is typically 
much later, unless the ISM density is very large and compa- 
rable to that of giant molecular clouds in the Galaxy (Dai 
& Lu 1999). 

The post-shock electrons and magnetic fields are usually 
assumed to share the fractions e e and es of the total internal 
energy, respectively. Therefore, the magnetic field in the co- 
moving frame is 

B' = y/32TTe B m p nyc w OA-ye^n 1 ' 2 G, (6) 

while the minimum Lorentz factor of accelerated electrons 
is (Huang et al. 2000) 

7m = e e ^^(7-l) + l, (7) 
p — 1 m e 

where m p and m e are the proton and electron masses, and p 
is the index of shock-accelerated electron distribution. The 
electrons can be cooled down within the dynamical time 
scale. The cooling Lorentz factor of electrons by purely syn- 
chrotron cooling is (Sari et al. 1998) 

syn 6Trrn e c(l + z) 
7c = a^B'H ' (8) 
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where c is the speed of light and ctt is the Thomson cross sec- 
tion. Further cooling of electrons by inverse-Compton (IC) 
up-scattering of the primary synchrotron photons can re- 
duce the above cooling Lorentz factor (Waxman 1997; Wei 
& Lu 1998; Panaitescu & Kumar 2000). Following Sari & 
Esin (2001), we get the final cooling Lorentz factor 



7c = 



1 + 4- 




(9) 



in the fast-cooling phase (7c < 7m), and 



„,syn 
£B 7c s,2/(4-p) 

; ) lm, 



<(-) 



e e s,l/(p-2) 



7c 



7m, 



7c > I — 



e Nl/(p-2) 



(10) 



in the slow-cooling phase (j m < j c ). In deriving the above 
7c we adopt 2 < p < 4, which seems to be easily sat- 
isfied by both the theoretical prediction and observations 
(Wu, Dai, & Liang 2004). It can be seen from equations 
(9) and (10) that y c always equals to ■yl vn for e e < es. In 
the case of e e > e B , the evaluation of j c is considered in 
three evolutionary stages. Equation (9) describes the first 
IC-dominated fast-cooling case, while equation (10) rep- 
resents the following IC-dominated and then synchrotron- 
dominated slow-cooling cases. The two characteristic fre- 
quencies of synchrotron radiation in the co-moving frame are 

eB' 2 



the typical frequency v'^ 
eB' 



2irm e c 



7 m and the cooling fre- 



quency v c = 7 C , which corresponds to the electrons 

27rm e c 

with minimum Lorentz factor y m and with cooling Lorentz 
factor 7c, respectively. Sari et al. (1998) showed that the co- 
moving peak spectral power per electron is independent of 



the electron energy, i.e. P' v i . 



B' , where e is the 



v .max 

electron charge. The co-moving peak intensity is therefore 
(Rossi et al. 2004) 



v .max 



K'.max nR m e c 2 a T „, _ 
-B nR. 



(11) 



4-7T 3 367re 
The spectral energy distribution of synchrotron radia- 
tion can be constructed by the characteristic frequencies v' c , 
v' m and the peak intensity I' v , max . In the fast-cooling phase 
(y' c < z/ m ), we have 



V , -V, x < 

^-u' J-v 1 .max ^ \ 



(77)" 



V < v c , 

v' c < v < v' m , (12) 

, v'm < V', 

and in the slow-cooling phase (y' m < u' c ), we have 

V' ,1/3 



l^m 



V , = V , x < 

± v t v ' ,max ^ » 



y v' ' 

V c s-(p-l)/2/f N-p/2 



V < V m , 

v' m < v < v' c , (13) 

, V c < V . 



We do not include the synchrotron-self-absorption effect 
which may be important mainly at the radio wavelength. 
Since afterglow polarizations had been detected mostly in 




Figure 1. Projection of a two-component jet. The inner cone 
is the narrow component, and the outer hollow cone is the wide 
component. The viewing angle between the jet axis (C.A.) and 
the line of sight (LOS) is 9y. Also shown are the apertures relative 
to the origin of the GRB jet, with ^ of the narrow component 
and 9w of the wide component. An arc centered to LOS with 
aperture 9 (dotted line) is divided into three segments, with the 
segment (xiv,27r — xjv) in the narrow component and the segments 
(Xw-iXn) an d (2?r — xat^tt — Xw) 111 the wide component. 



optical bands, we will only investigate the temporal evolu- 
tion of optical polarization. 



3 POLARIZATION 
3.1 Formulation 

The polarization of emission from a point-like region aver- 
aged by the distribution of the post-shock tangled magnetic 
field is (Sari 1999; Gruzinov 1999) 



11(7,0) = n 5 (7, 0) sin e 



e-i 



2e + (l-C 2 )5 2 (7,e)sin 2 0' 



(14) 



where 6 is the angle between the velocity vector of the emit- 
ting point and the line of sight, 5(7, 6) = 



1 



7(1-/3 cos (9) 

is the Doppler boosting factor. The parameter £ 2 = 
(B^) /2{B 2 ) denotes the level of anisotropy of the magnetic 
field distribution, where B± and B\\ are magnetic field com- 
ponents perpendicular and parallel to the normal of the 
shock plane. The triangular brackets mean that the quan- 
tity inside is averaged over the solid angles, no is the lin- 
ear polarization degree of synchrotron photons emitted by 
the electrons in an ordered magnetic field. For an isotropic 
power- law distribution of electrons with index p, the approx- 
imate no = (p + l)/(p + 7/3) holds true for a wide range of 
frequencies. The exact dependence of Ilo on the frequency 
is given in the Appendix (see Figure Al). Using the above 
approximate expression for no actually leads to neglectable 
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errors in the final results. We therefore adopt the usual ap- 
proximate value of IIo in the following calculations. 

To determine the total polarization of a GRB jet, we 
need to integrate the Stokes parameters over the jet surface. 
From Figure 1, by virtue of spherical geometry we have 



Xi(0) 



■ — cos 1 ( 



ire{0 v -0j) 

cos 0j — cos 8 cos 8v 



sin ( 

Ti 



sin V 



,0<0-, 
0- < 
0+ < 



<0+, (15) 



where 0y is the observer's viewing angle with respect to the 



jet axis, 0- 



and #4 



+ 



Q(x) is the 



heaviside step function with Q(x > 0) = 1 and Q(x < 0) = 

0. Note that equation (15) is valid for both components, 

1. e. Xi — XN when 0j — 8n, and Xi = Xw when 0j = 
8w ■ Ghisellini & Lazzati (1999) have given an approximate 
expression of Xi (ipi in their work) for the intermediate case. 
Integrating over each component of the jet, we get the Stokes 
parameters as 



r 9 v +9 w 



Qi(t 



^) = ^2 £ / S 3 (^,0)ll,M^,0)RUin8d8 

L Jo 

rxi, max 

x / cos2ipd(p, (16) 



1 4- 7 f 0v+f>w 

Ii{t,v) = -^- / 5 3 {~ fi ,0)l' v ,, i R t i s™0d0 



L JO 



- r 



dtp, 



(17) 



where i — W or N, v' — (1 + z)5 1 {'fi,0)v is the fre- 
quency in the co-moving frame. For the narrow component, 
XN,min = Xn(0) and XN,max = 2n - xn(8)- For the wide 
component, the azimuthal integral is from xw,min = Xw{0) 
to xw,max = Xn,o(0), and again from xw,min = 27r-Xiv,o(#) 
to xw,max = 27r — xw (8) ■ Here Xn,o(0) is calculated by 
inserting the inner half-opening angle of the wide hollow 
cone 0n,o into equation (15), where 0n,o is the initial aper- 
ture of the narrow component. Note that Stokes parameter 
Un = Uw = can be easily derived from the symmetry. The 
Stokes parameters for the whole jet are therefore calculated 
by summation of the values of these two components, 



Qtotit, v) 



l + Z 



n 2 , 

L Jo 



V +6w 



smddd x 



{ 8 3 (-yw ,8)ll' ,w^l(lw ,8)Rw(sin2xN,o — s\n2xw) 
-S 3 ( lN , 0)l'„, N Il{ lNl 0)R 2 N sin 2 X n}, (18) 



Qtot(t,u) 



l + Z 

D 2 , 



sm0d8 x 



{ 5 3 (7i4/,6 l )/^ iH/ n(7i V ,6 l )i?^(sin2xjv,o - s\n2xw) 
-5 3 (-y N , 6>)7^, Ar n( 7J v, 0)R 2 N sin 2 X at}, (19) 

Itot{t,v) = —J- / {5 3 ('yw,8)l' l ,>,wRw('2XN,o-2xw) 



D 2 \ 

L Jo 

+S 3 ( 7JV , 0)l' v , tN R 2 N (2tt - 2 X iv) } sin 0d0. (20) 

The total polarization of the two component jet is thus given 
by 



Ptot — 



Qt 



(21) 



The corresponding polarization vector is either (anti-) par- 
allel or perpendicular to the vector from C.A. to LOS pro- 
jected in the sky (see Figure 1), depending on whether the 
magnetic field parameter £ 2 is larger than unity or not, as 
well as on which direction the observer views and how the 
jet evolves (Ghisellini & Lazzati 1999; Sari 1999). 



3.2 Numerical Results 

Below we use several sets of physical parameters to illustrate 
the temporal evolution of afterglow polarization of a two- 
component jet. The main interesting parameter effects may 
arise from the ratio of the isotropic-equivalent energy of the 
two components Ew,iso/ En,xso, and the ratio of their initial 
half-opening angles 8w,o/0n,o- We fix other parameters and 
assume En,iso = 10 53 erg, 0n,o = 0.05, n = 1 cm~ 3 , e e = 
0.1, <lb = 0.01, p = 2.2, £ 2 = 9 and n = 60% in the 
following calculations. The equal arrival time surface effect is 
neglected in this primary research. The adopted parameters 
of the cosmology are Q m = 0.27, Q.a = 0.73, Ho = 71 km s _1 
Mpc -1 , and the GRB is assumed to be located at z — 1.0. 

Figure 2 shows the light curves and polarization evo- 
lutions of GRB optical afterglows (u opt = 4.5 x 10 14 Hz) 
from two-component jets 2 , which are assumed to evolve 
with no lateral expansions. The ratio of the wide component 
aperture to that of the narrow component is chosen to be 
0w,o/0n,o = 6. The upper panel corresponds to a relatively 
low isotropic energy ratio as Ew,iso/ 'En ,iso = 0.1, while the 
lower panel corresponds to a high ratio as Ew,iso/ En,iso = 
0.3. When the line of sight is in the narrow component, i.e. 
0v < 0n,o, the light curve differs little between each other 
since the light mainly comes from the narrow component. 
The light curve steepens from ~ t -0 9 to t~ 16 at around 
tj,N ~ 0.2 day, which is shallow as expected for the non- 
lateral expansion jet (Meszaros & Rees 1999). The polariza- 
tion degree is also dominated by the narrow component, and 
changes from an initially negative value to a later positive 
value, which means that the position angle of polarization 
changes by 90°. After reaching a maximum, the polariza- 
tion decreases and approaches zero. When the line of sight 
lies within the wide component, the light curve flattens when 
the narrow component decelerates and enters into the field of 
view. The degree of flattening is determined by the viewing 
angle 0v and the energy ratio Ew,iso/ 'En tiso . The smaller 
0v or Ew,iso/ EN,iso is, the more contribution is resulted 
from the narrow component, as can be seen in the figure. 
Contrary to the light curve, the total polarization follows 
almost the evolution of polarization of the narrow compo- 
nent for a long time. This is because the narrow component 
is off-viewed and the increase of the point polarization is 
greater than the decrease of the Doppler boosting when 
increases (equations 14 and 16). The only exception is when 
0v/8n,o = 5, in which case the polarization evolution is 
dominated by the wide component. For the other cases, the 



2 To compare which component dominates the light curve or po- 
larization, we treat the two components separately, and then add 
them together. In the figures we only give the total flux density 
and polarization. 
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Figure 2. Light curves and polarization evolutions of GRB op- 
tical afterglows from two-component jets. The jet is assumed to 
have no lateral expansion. The contrasts between the wide and 
narrow component parameters are Eyy i so / E jv i so = 0.1 (upper 
panel) or 0.3 (lower panel), and 8w,o/^N,o = 6. Different lines 
correspond to different viewing angles (By). 

polarization level is depleted to a certain extent by the wide 
component. In the cases of 9v/9n,o = 1.2 and 1.5, there 
exist long periods with nearly constant polarization levels 
and unchanged position angles, which is the typical char- 
acteristics of off-viewed jets (Granot et al. 2002; Rossi et 
al. 2004). In general, increasing the contrast between Ew,iso 
and En,xso would decrease the maximum of polarization de- 
gree and reduce the flattening or bump by the central narrow 
component. 

Figures 3 and 4 depict the light curves and polarization 
evolutions of GRB optical afterglows from two-component 
jets with lateral expansions. When 9y < 9n,o, the light curve 
is dominated by the narrow component and experience a 
steepening from t~ ' 9 to t~ 2 ' 2 at tj t tf « 0.2 day. The light 
curve flattens at late times when the wide component be- 
gins to dominate over the narrow one for 6w,o/Sn,o = 6, 
in which case the jet break of the wide component occurs 



Figure 3. Light curves and polarization evolutions of GRB op- 
tical afterglows from two-component jets. The jet is assumed to 
expand laterally at the speed of light. The contrasts between the 
wide and narrow component parameters are E\yi so / Ejy ; so = 0.1 
(upper panel) or 0.3 (lower panel), and 0\y,o/6N,a = 2. Different 
lines correspond to different viewing angles (9y). 

at rather late times, tj,w > 10 days 3 . The sign of polar- 
ization changes twice for 8y < &n,o- In most cases, except 
when 8y is very close to 0\y,o, the polarization curves have 
one peak at around ij,jv. These are typical for jets with lat- 
eral expansion (Sari 1999). For 8y > 9n,o, the polarization 
is mainly dominated by the narrow component before ij.jv- 
Exceptions are 9y — 1.88n,o in the case of 9w,o = 29n,o and 
9y — 59n,o in the case of 9w,o = 69n,o, where the polariza- 
tion is dominated by the wide component during the whole 
evolution. After the jet break the asymmetry of the narrow 
component diminishes due to sideways expansion, the po- 
larization begins to be determined by the wide component, 
almost regardless of whether 9y is larger than 9n,o or not. 

3 There is another possibility to cause a late time flattening, in 
which case the jet (single component) expanded sideways and 
becomes spherical and non-relativistic (Huang & Cheng 2003). 
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Figure 4. Same as Figure 3 except for 9w,o/^N,o = 6 now. 

In the figures we can see that the absolute polarization level 
at t > tj,N increases with Ew,iso/ En,xso- As shown in the 
figures, the light curves in cases of 9y > 0n,o are influenced 
by the narrow component through a bump or flattening. 
When the ratio 8w,o/8n,o or especially Ew,iso/ En,iso in- 
creases, the influence of the narrow component to the light 
curve is reduced, and the maximum polarization level also 
decreases. 

There are mainly two differences between the scenario 
without lateral expansion and that with lateral expansion. 
The first is the temporal behavior of light curve after the 
jet break time tj . The light-curve steepening of a non-lateral 
expansion jet around tj is relatively shallow, while the steep- 
ening of a lateral expansion jet is large. The second is the 
width of the peak in the temporal polarization curve. A 
jet with lateral expansion has a more narrow polarization 
peak than a jet without lateral expansion, as can be seen by 
comparing Figure 2 with Figure 3. Rossi et al. (2004) have 
investigated the polarization evolutions of structured jets, 
with both the power-law and Gaussian distributions. They 
found that the position angle of polarization is not changed 
in time for these structured jets. A two-component with- 



out lateral expansion also has an unchanged position angle 
when 9m,o < @v < &w,o ■ The similarities between such 
two-component jets and those with power law or Gaussian 
structure are that they all have no or little lateral expansions 
and they are more energetic near the jet axis. 

3.3 The case of GRB 020813 

The bright long duration GRB 020813 detected by HETE- 
II spacecraft had been well-observed for evolutions of both 
optical flux density and polarization. The redshift of this 
burst is z = 1.255, by identifying the [O II] A3727 line of 
the host galaxy (Barth et al. 2003). The light curve of GRB 
020813 afterglow is the smoothest one among GRBs local- 
ized so far (Laursen & Stanek 2003). The optical afterglow 
has a jet break around 0.33 - 0.88 day since the main burst 
(Gorosabel et al. 2004; see also Covino et al. 2003b, Li et al. 
2003). The temporal index before the break (ai) is ~ —0.55 
— ~ —0.78, and the index after the break (02) is ~ —1.44 

1.75. Barth et al. (2003) have obtained the afterglow 

spectrum in optical band, and found that the spectral index 
is f3 = 1.06 ± 0.01. Assuming that the spectrum is fast cool- 
ing, the inferred power law index of electron distribution is 
p = 2.12 ± 0.02. The spectropolarimetric observations had 
been made by Keck between 0.19 — 0.33 days and by VLT 
between 0.88 - 4.05 days (Barth et al. 2003; Gorosabel et 
al. 2004). The degree of polarization decreases from about 
3% to less than 1%, while the position angle varies slowly 
and can be regarded as constant. Both observations showed 
that the majority of these detected polarizations are intrin- 
sic and not effected by the line-of-sight dust. Therefore, the 
afterglow of GRB 020813 provides us an ideal sample to 
investigate the structure of GRB jets (Lazzati et al. 2004). 

We have fitted the V-band light curve and polarization 
evolution of GRB 020813 afterglow in the context of the 
two-component jet model. The data of light curve and po- 
larizations are taken from Gorosabel et al. (2004), in which 
they included and averaged the polarization observations 
made by Barth et al. (2003). We assume that the position 
angle of polarization is unchanged, at least within the time 
scale we concerned. The lateral spreading jet scenario can 
be excluded for GRB 020813 because the observed temporal 
index Q2 is larger than the inferred —p. Even adopting a 
harder electron distribution with 1 < p < 2, the predicted 
0L2 = -(p + 6)/4 is smaller than -1.75 (Dai & Cheng 2001). 
However, a non-lateral expansion jet running into the ISM 
is able to reproduce the observed afterglow light curve. In 
this case the theoretical Aa = ai — a.2 = 0.75 is consistent 
with the observed value, ranging from Aa b s ,v = 0.50±0.17 
(Covino et al. 2003b) to Aq o6s ,v = 1.20 ± 0.17 (Gorosabel 
et al. 2004). The less constrained Aa„i, Si v is due to the less 
constrained jet break sharpness. Since the detected position 
angle of polarization did not change around the jet break 
time ttreak,v = 0.33 — 0.88 days, we expect the line of sight 
is located within the wide component (e.g., see Figure 2 for 
the case of By = 1.2#/v,o). However, by analyzing the polar- 
ization evolution, the viewing angle should be close to the 
narrow component. Thus the wide component is responsi- 
ble for the main burst and early afterglow, while the nar- 
row component is for the whole evolution of the polarization 
and the late (mid) time afterglow. The fittings to the GRB 
020813 V-band afterglow including the light curve and po- 
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Figure 5. Fitting to the GRB 020813 V-band afterglow light 
curve (upper panel) and its polarization evolution (lower panel). 
The Galactic extinction in the V-band Ay = 0.42 is adopted 
(Covino et al. 2003b). 

larization evolution are shown in Figure 5, where parame- 
ters of the two-component jet are Eu,iso = 1.2 x 10 54 erg, 
9m = 0.02, Ew.iso = 0.2i?jv,iso, and 6w = 3.50n- The jet is 
viewed at 6v = 1.450jv. The shock parameters are e e = 0.05, 
e B = 5 x 10 -5 , p = 2.15 and £ 2 = 1.25. The ISM density is 
n = 1.0 cm -3 . The relatively low £ 2 was also indicated by 
Barth et al. (2003). Granot & Konigl (2003) proposed that 
the magnetic field configuration may be isotropized by the 
post-shock turbulence, and £ 2 we observed is the averaged 
value over the shell width. 

Due to the sparse detections of polarizations, especially 
at t ~ 0.5 day, Barth et al. (2003) argued that it is still 
possible to accommodate the lateral spreading jet model for 
GRB 020813. But a structured jet, with the center to be 
brighter or more energetic while the wings to be dimmer or 
less energetic, seems more natural as compared with other 
models according to polarization observations (Lazzati et 
al. 2004). However, we should also bear in mind that the 
details of the structures, i.e., whether its profile is power- 
law, Gaussian, or two-component, still cannot be definitely 
determined currently. 



4 CONCLUSIONS AND DISCUSSION 

In this paper we have investigated the polarization prop- 
erty of a two-component GRB jet. The inner component is 
narrow and more energetic while the outer one is wide and 
less energetic. The hydrodynamic evolution of each compo- 
nent follows the relativistic Blandford-McKee scaling law in 
the ISM case. The shock physics is assumed to be the same 
for these two components. We include the inverse Compton 
scattering in the calculations of afterglow light curves. The 
major effect of inverse Compton scattering is to further re- 
duce the cooling Lorentz factor of electrons. Since the typical 
energy of inverse Compton scattered photons lies above the 
X-ray band, and the synchrotron-self absorption operates 
in the radio afterglow, we considered only the synchrotron 
origin of optical polarization. 

The resulting optical light curves and polarization evo- 
lutions depend strongly on the ratios of the intrinsic pa- 
rameters of the wide component to those of the narrow one 
(i.e. Ew.iso/ E N , mo and 9w,o/Qn,o; see also Peng, Konigl, & 
Granot 2004, in which they investigated the light curves 
of GRB afterglows within the two-component jet model in 
analytical way), on the different assumptions of the lateral 
expansion, and on the observer's viewing angle as well. Two 
scenarios of different lateral expansions are studied. One is 
to assume that the jet has no lateral expansion, while an- 
other is to assume the jet expands laterally at the speed of 
light. In both cases we find that within a wide range of the 
viewing angle the polarization is dominated by the narrow 
component for a long period. For an observer with the line 
of sight located in the wide component, the position angle 
of polarization will not change if the jet has no lateral ex- 
pansion. This is because the geometric asymmetry of the 
narrow component does not change with respect to the line 
of sight. On the other hand, the position angle will change 
by 90 degree for a laterally expanding jet since the jet ge- 
ometry evolves from an asymmetrical one to a symmetrical 
one, with the transition taking place around tj (Ghisellini & 
Lazzati 1999; Sari 1999). Although the two-component jet 
model has been proposed to explain some peculiar afterglow 
light curves (Frail et al. 2000; Berger et al. 2003; Huang et 
al. 2004), it does not imply necessarily that the light curve 
of a two-component jet must be peculiar. When the line of 
sight lies within one component, another component may 
cause a flattening of the light curve at late times, not neces- 
sarily a bump, depending on the contrast of the parameters 
between these two components (see also Huang et al. 2004, 
where the equal arrival time surface effect was included). 

We find that the two-component jet model with the 
assumption of no lateral expansion is able to explain the af- 
terglow of GRB 020813 through fittings to its optical light 
curve and polarization evolution. The essence of this ex- 
planation is to assume the line of sight is located in the 
wide component, which ensures a constant position angle of 
polarization. However, as pointed by Lazzati et al. (2004), 
any jet with a more energetic core and less energetic wings 
will produce qualitatively similar polarization curves. Due 
to the sparse polarization data of GRB 020813 afterglow, 
especially at very early times and around the jet break time 
tj ~ 0.5 day, we cannot determine the actual structure of 
GRB jets definitely. In the upcoming Swift era, we expect 
that several well-observed GRB afterglows with copious data 
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of both magnitude and polarizations will give us new clues 
about jet structures. 
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APPENDIX A: EVALUATION OF THE 
SYNCHROTRON n IN AN ORDERED 
MAGNETIC FIELD 

The linear polarization of optically thin synchrotron radia- 
tion in an ordered magnetic field is calculated by averaging 
over an isotropic distribution of electrons (Longair 1994) 



IlofV) 



J G(x)N(-y e )d>y e 
/ F{x)N{ le )d le ' 



(Al) 



where x 



with v' lt 



47vm e c ' 

x J K 5 / 3 (t)dt and G(x) = xK 2 /3(x), where Ks/3(x) and 
K2/3 are the modified Bessel functions. In the limit of x <C 1, 

F(x) » 2G(x) « ^ . (^) 1/3 , while in the limit of 



%/3r(l/3) 2 



x > 1, F(x) « G(x) 



The electron distribu- 



ter 1/2 - 
-x' e 

tion function, N(-y e ) = — — , determines the dependence of 

IIo on the frequency v' . 

Fast-cooling case (<y c < 7m )- The electron distribution is 
approximated by a broken power-law, with N(-y e ) oc 7 < T 2 for 
7c < 7e < 7m an d N{j e ) oc 7 ( r p ~ 1 for -y m < 7e < 7M , where 
7 m is the maximum Lorentz factor of accelerated electrons. 
Rewriting equation (Al), we get 



n (i/) 



\H G^X^dx + X^^ HZ G{x)x { V-^l 2 dx 

£ C F(x)x~y 2 dx + x£- p)/2 f* m F{x)x(p- 2 )/ 2 dx ' 



(A2) 



where x c = u /v' c , x m = v' jv' m and xm = v /u M , Flo ap- 
proaches 1/2 for v' <C v' c f' m , 9/13 for v' c ^ v' <C v' m , 
and (p + 2)/(p + 10/3) for u' m < u' < v' M . 

Slow-cooling case (-y m < 7c ). The electron distribution 
is approximated by a broken power-law, with iV( 7e ) oc 7,T P 
for 7m < 7 e < 7c and iV(7 e ) oc ^ p ^ 1 for 7c < Is < 1m- 
Equation (Al) can be deduced to be 



W) = 



£ c m G(x)x^ 3 ^ 2 dx + x- 1/2 f* C M G{x)x {p - 2 ^ 2 dx 
\ Xm F(x)x(p- :i )/ 2 dx + Xc 1/2 f Xc F(x)x(p- 2 )/ 2 dx 



.(A3) 



which approaches 1/2 for v' <C v' m <<iC v' c , (p + l)/(p + 7/3) 
for v' m < v' < v' c , and (p + 2)/(p + 10/3) for v' c < < 1/^. 

The above approaches are consistent with those in Ta- 
ble 1 of Granot (2003), in which the simple analytic expres- 
sion for flo as a function of the spectral slope for optically 
thin synchrotron emission can be found. Figure Al clearly 




Figure Al. The dependence of IIo on the observed frequency v. 
Panel (a): the fast-cooling case {y c < u m )', Panel (b): the slow- 
cooling case (y c > u m ). The dashed, solid and dotted lines corre- 
spond to p = 2.1, 2.3, and 2.5 respectively. 



shows how the polarization degree Bio varies with the ob- 
served frequency v — S(j, 8)(l + z)~ x v' . As can be seen, in 
both fast and slow cooling cases, there are three platforms 
in each polarization curve, which represent the three typical 
polarization levels discussed above. 
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